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ABSTRACT 

The ionizing ultraviolet background (UVB) during reionization can suppress the gas content 
of low-mass galaxies, even those capable of efficient atomic cooling (i.e. with virial temper- 
atures T V ii > 10 4 K). This negative radiative feedback mechanism can thus reduce the star 
formation efficiencies of these halos, which may delay the completion of reionization. In this 
work, we explore the importance of UV radiative feedback on T V j r > 10 4 K halos during the 
middle and late stages of reionization. We do not try to self-consistently model reionization; 
instead, we explore a large parameter space in an attempt to draw general, robust conclusions. 
We use a tiered approach. Using 1-D hydrodynamical simulations, we model the ability of gas 
to collapse onto halos of various masses under UVBs of various intensities. We then gener- 
ate realistic, parametrized maps of the inhomogeneous UVB, using large-scale semi-numeric 
simulations. By combining these results, we find that under all reasonably conservative sce- 
narios, UV feedback on atomic ally-cooled halos is not strong enough to notably delay the 
bulk of reionization. Such a delay is only likely if ionizing efficiencies of z > 10 sources are 
much higher (~ two orders of magnitude) than z ~ 6 data seem to imply. Towards the end 
of reionization, star formation can be quenched only in halos in a narrow mass range close to 
the atomic-cooling threshold. This result depends only weakly on the intensity of the UVB: 
quenching star formation in halos just twice as massive requires an order of magnitude in- 
crease in source ionizing efficiencies. This implies that the natural time-scale for the bulk of 
reionization is the growth of the global collapsed fraction contained in T vlr > 10 4 K halos. 
Thus the likely reionization scenario would involve a small HII filling factor "tail" extending 
to high redshifts, governed by more complicated feedback on T vu < 10 4 K objects, followed 
by a period of relatively rapid evolution in the HII filling factor. Furthermore, our results 
underscore the importance of extended dynamical ranges when modeling reionization. Sim- 
ulations must be capable of resolving halos with mass > 10 8 M Q , even when modeling the 
late stages of reionization, while at the same time being large enough to capture HII regions 
several tens of Mpc in size. 

Key words: cosmology: theory - early Universe - galaxies: formation - high-redshift - 
evolution 



1 INTRODUCTION 

The epoch of reionization, when light from early generations of as- 
trophysical objects began flowing through the intergalactic medium 
(IGM), remains one of the most compelling frontiers of mod- 
ern cosmology, offering a wealth of information about cosmo- 
logical structure formation and physical processes in the early 
universe. Only recently have we begun to gather clues concern- 
ing this epoch. The presence of flux in the Lyman line absorp- 
tion regions of z < 6 quasars discovered in the Sloan Digi- 
tal Sky Survey (SDSS) indicates that reionization is mostly com- 
pleted by this redshift (though note that large-scale fluctuations 
in the ionization field can yield long ionized sightlines even be- 
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fore the completion of reionization, e.g. ILidz et al.|[2007h . Spec- 
tra from some higher redshift quasars have even been inter preted 
as evidence of a partially neutral IGM (Mesinger & Haima ril 20041 : 
IWvithe & Loebll2004l: iFan et al.|[200r3 ; iMesinger & Haimanll20070. 
though most of these interpreta t ions are controversial dLidz et al. 
1 20061 : [Bolton & Haehneltl |2007|; iMaselli et al.l [20071 : iBecker et~ 
120071 ; IMesinger & Furlanet to 2008a). The Wilkinson Microwave 
Anisotropy Probe (WMAP) meas ured a Thomson sca ttering opti- 
cal depth of r e = 0.087 ± 0.017 jPunklev et ai]|2008l) . Assuming 
instantaneous reionization, this value of r e corresponds to a reion- 
ization redshift of z — 11.0 ± 1.4. However, this is only an inte- 
grated measurement and can tell us very little about the evolution 
of reionization. 

With the interpretation of even the sparse existing data being 
subtle and complicated, much effort has been invested in improving 
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the modeling of reionization. We think we can accurately model the 
abundances and clustering properties of dark matter halos at high 
redshift. Combined with some analytic relation of halo mass — > ef- 
fective ionizing photon emissivity, one can get an estimate of how 
the global neutral fraction evolves; this technique has been used 
is countless reionization studies. Aside from the seemingly magi- 
cal "— »" in the previous sentence, one of the main additional un- 
certainties are the feedback processes: how do sources impact the 
current and future generations of sources? Answering this question 
is non-trivial, especially in the early universe. Radiative, chemi- 
cal and mechanical feedback can effect the ionizing efficiencies 
of the first, highly biased sources. Molecular hydrogen (H2) cool- 
ing can allow very small (with virial temperatures, T v i r , of several 
hundred Kelvin) halos to host astrophysical s ources. At later time s 
(though still prior to the bulk of reionization, lHaiman et al.ll 19971) . 
as the H2 dissociative background builds-up and the contribution 
of molecularly-cooled halos declines, feedback enters the regime 
of atomically-cooled halos, T v i r > 10 4 K. This regime is more 
straightforward to model as negative radiative feedback emerges as 
a single dominant mechanism, especially on large scales. 

In the presence of an ionizing background radiation, the 
IGM is photo-heated to a temperature of > 10 4 K, raising 
the cosmological Jeans mass, w hich could suppr e ss gas accre- 
tion onto small-mass halos (e . g. lEfstathioul 1 19921: Ishapiro et al.l 



1 1994 iThoul & Weinberg! IT996; Hui & Gnedmlll997l) . Early work 
on this subject (so-called "Jeans mass filtering") suggested 
that an ionizing background would completely suppress star 
formation in low-redshift "dwarf galaxy" halos with circu- 
lar velocities w C i rc < 35 km s" 1 , and partially suppress 
star-formation in halos with 35 km s" 1 < v c i rc < 100 
km s~ - 1 jThoul & Weinberdll996l) . Many reionization studies as- 
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IWvifhe & Cenll200l llliev et alj|2007l; IWvithe et al.ll2008l). How 



ever, more recent studies jKitavama & Ikeuc hi 2000; Dijkstra et al.l 
|2004|) find that at z > 3, more compact halo profiles, increased 
cooling efficiencies, and shorter exposure times to the ultraviolet 
background (UVB) could lessen the importance of negative radia- 
tive feedback. We expand on these works by exploring a wider 
parameter space, placed in a broader context of an inhomoge- 
nious reionization with a patchy UVB du e to both source clustering 
and modulation by HII regions. iGnedirj J2000h studied the effects 
of such radiative feedback in cosmological simulations; however, 
these simulation boxes were of necessity very small and only in- 
cluded a single reionization model. 

The purpose of this paper is to explore the importance of 
UV radiative feedback during the middle and late stages of reion- 
ization. We do not attempt to self-consistently model feedback 
during reionization; such a thing is beyond the capacity of cur- 
rent simulations, especially given our poor understanding of the 
first generations of astrophysical sources and their environments. 
Instead we statistically present the effects of an inhomogeneous 
UVB on the suppression of gas content in low-mass galaxies ca- 
pable of atomic cooling, i.e. T v i r > 10 4 K. We do this using a 
tiered approach: using nu merical simulations jThoul & Weinberg 
1 1 9951 : |Pii kstra et al.l2004l) to calibra te very large scale, high resolu - 
tion "semi-numerical" simulations dMesinger & Furlanettoll2007l) . 
In order to keep things general, we explore wide swaths of param- 
eter space keeping assumptions minimal. 

This paper is organized as follows. In ^2. 1 1 and i]2.2| we de- 
scribe our hydrodynamic collapse simulations and semi-numerical 



cosmological simulations, respectively. In £13.11 we present the re- 
sults from our collapse simulations covering a wide range of pa- 
rameter choices, while in i)3.2| we present parametrized distribu- 
tions of UV fluxes from our cosmological simulations. In S]3.3l we 
combine these results to quantify the importance of UV feedback 
during reionization. In fj4]we discuss the assumptions and uncer- 
tainties in our approach. Finally, in Sj5]we present our conclusions. 

We quote all quantities in comoving units, with the exception 
of flux, where we denote proper units with a prefix 'p'. We adopt 
the background cosmological parameters (J1a, ^m, ^6, n, og, Ho) 
= (0.76, 0.24, 0.0407, 0.96, 0.76, 72 km s - 1 Mpc" 1 ), matchin g 
the three-year results of the WMAP satellite dSpergel et ai]l2007l) . 
which in turn are con sistent with the recent five-year data release 
dKomatsu et al]2008h . 



2 SIMULATIONS 

As we outlined above, we use spherically-symmetric simulations 
to model the collapse of gas onto halos with masses M > 10 8 Mq . 
With these numerical simulations, we are able to determine the 
amount of gas which collapses onto a halo of mass M at z, un- 
der the presence of a UVB with specific intensity J21 turned on 
at redshift z on . We then generate large-scale halo, ionization and 
parametrized flux fields, using semi-numerical techniques. Com- 
bining these results, we investigate what fraction of halos experi- 
ence strong negative radiative feedback. We describe our method- 
ology and numerical techniques in greater detail below. 



2.1 Hydrodynamical Collapse Simulations 

We investigate the impact of a photoionizing flux on the abil- 
ity of gas to cool and collapse onto dark matter halos by per- 
forming 1-D hydrodynamical simulations which include both gas 
and dark matter. We use the code that was originally written by 
IThoul & Weinberg] dl995h and modifiedcharacterized to study the 
impact of photoionization feedback on the formation of high- 
redshift low-mass galaxies (see Dijkstra et al 2004 for a more de- 
tailed description). 

To summarize our calculations: an externally generated UVB 
impinges upon a spherically symmetric collapsing dark matter halo 
of mass M that would have collapsed to r — at redshift z, in the 
absence of gas pressure. The ionizing radiation field is switched on 
at z on , and is characterized by a power-law spectrum of the form 
J{v) = Ji\{vlvHY a x 10~ 21 ergs" 1 Hz" 1 pern" 2 sr" 1 , where 
vh is the ionization threshold of hydrogen. We are interested in the 
fraction of the gas, / gas = M gas (p)/Af gas (p = 0), (i.e. the mass 
of collapsed gas divided by the mass of gas which would have col- 
lapsed in the absence of pressure) that is instantly available for star 
formation at z. Therefore, we compute the fraction of gas which has 
collapsed at 2 (with an associated Hubble time t C oii) which differs 
from the calculations presented by e.g. Thoul & Weinberg (1996), 
and Dijkstra et al. (2004) who computed the fraction of gas which 
was able to collapse up to 2t co n. 

For all runs, we switch on the UVB at z on = 14. The 5- 
yr WMAP polarization data place a constraint on the reionization 
redshift (assuming instantaneous reionization) of z = 11. Oil. 4 
dDunklev et al . 2008). Thus it is highly unlikely that the majority of 
the universe was ionized before z ~ 14 (for so me interpretations 
of the less-constrained 3-yr data, s ee, for example lHaiman & Brvanl 
2006; IChoudhury & Fe rrara 2006). Therefore, our choice of z on — 
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14 is conservative, since most halos should be exposed to ioniz- 
ing radiation at later times and so we overestimate the impact of 
radiative feedback. 

Furthermore, following Dijkstra et al. (2004) we have as- 
sumed a = 1.0, which is typ ically associated with accreting black 
holes (e.g. Madau et al. 2004) and not with st ellar sources which 
would be better characteriz ed with a ~ 5 dThoul & Weinberg! 
ll996l : lBarkana" & Loebll200lh . Because the simulations ignore self- 
shielding, lowering a is equivalent to increasing J21 in terms of the 
photoionization and photoheating rates. Therefore, in our calcula- 
tions we overestimate the impact of the photoionization feedback 
(see §|4]for a more extended discussion). 



2.2 Semi-Numerical Cosmological Simulations 

We use an excursion-set approach combined with first-order La- 
grangian perturbation theory to efficiently generate density, veloc- 
ity, halo, and ionization fields at z =7, 10 and 13. This "semi- 
numer ical" simulation is presented in Mes inger & Furlanettd 
d2007l) . to which we refer the reader for details. A similar halo- 
finding scheme has also been presented bv lBond&Mversl ( ll996l) 
and a si milar scheme to ge nerate ionization fields has been pre- 
sented bv lZahn et ail d2007h . 

Our simulation box is 100 Mpc on a side, with the final ion- 
ization fields having grid cell sizes of 0.5 Mpc. Halos with a total 
mass M > 2 x 10 7 Mq are filtered out of the linear density field 
using excursion-set theory; however, for this work, we only concern 
ourselves with halos capable of efficient atomic hydrogen cooling, 
T vir > 10 4 K, or M > 1O 8 M for the redshifts in question. Halo 
locations are then adjusted using first-order Lagrangian perturba- 
tion theory and velocity resolution of ~ 0.14 Mpc. The resulting 
halo field matches both the mass function a nd statistical clustering 
prope rties of halos in N-body simulations (Mes inger & Furl anetto 
I2OO7I) . 

In constructing the ionization field, the IGM is modeled as 
a two-phase medium, comprised of fully ionized and fully neu- 
tral regions (this is a fairly accurate assumption at high redshifts 
preceding the end of reionization, unless the X-ray background is 
rather strong). For each halo and density field at our three redshifts 
of interest, z = 7, 10, 13, we generate several ionization fields 
corresponding to different values of Shi by varying a single effi- 
ciency parameter, (, again using the excursion-set approach (c.f. 
iFurlanetto et al.ll2004h . 

This fast, semi-numeric approach is thus ideally suited to 
reionization studies, because we are able to "resolve" relatively 
small halos and simultaneously sample a large, representative vol- 
ume of ionized bubbles. 



2.2.1 Generating the Flux Field 

For each halo and ionization map, we create a corresponding UV 
flux field on a 200 3 grid, by summing the contributions of halos 
whose sightlines to the cell of interest go through ionized IGM. 
Specifically, we compute the flux of ionizing photons (in units of 
ionizing photons s _1 pcm -2 ) with 



f(x , 2) = (!+£)! y Mi x 6i °- e -i«-»»i/wp , 

47r ^r 1 x - Xi 2 



(1) 



where x is the location of the cell of interest, Mi is the total halo 
mass, Xi is its location, l m { p is the assumed mean free path (m.f.p.) 



of ionizing photons in the ionized IGM, and the factor of (1 + 2) 2 
converts the factor |x — Xi| 2 from comoving into proper units. We 
work in the Euclidian limit (i.e. the relevant distances are |x— xi | < 
m.f.p. <C hubble length, and cosmological corrections can be ig- 
nored; see eq.|2j. We assume Z m f p = 20 Mpc throughout, with one 
exception that is discussed explicitly in § 13,21 This value is consis- 
tent at z ~ 7 with extrapolations from the observed z < 4 abun- 
dances of Lyman limit systemJ](LLSs; e.g. Storrie-Lombardi et al.l 



ll994LlMiralda-Escudell2003l;IChoudhurv et alj2008h . Our choice is 
also consistent with the theoretical assumption that the m.f.p. in 
our ionized regions is dominated by self-shielded minihalos (halos 
with virial temperatures T v i r < 10 4 K). Assuming a roughly ran- 
dom distribution of minihalos wi th number density n m h and radius 
Rmh , this yields a m.f.p. of (c.f. lFurlanetto & Oh 2005): 



lmfp(z) 



7T n mh i? 2 lh 



17 ( M mh y f O05^ 



VlO 6 M 



0.12 



Mpc, 
(2) 

where Af m h is the typical minihalo mass, A m h its mean overden- 
sity, and / m h the fraction of matter locked up in minihalos. As de- 
fined above, / m h, can vary from 0.02 at 2 = 13 to 0.09 at z — 7. 
However, since the derivation of Z m f p and its evol ution (existing 
minih alos can get photoevaporated with time, e.g. IShapiro et al.l 
120041 partially countering the buildup of new minihalos) is ex- 
tremely uncertain, we assume a redshift-independent l m f p . We will 
elaborate more on how i m f p affects our conclusions below. 

Finally, and most importantly, ei on in eq. (Q3 denotes the rate 
at which ionizing photons are released into the IGM by a dark mat- 
ter halo per unit mass. This number depends on several factors: (i) 
the number of ionizing photons that is emitted per stellar mass (in 
Mq), e 7 , (ii) the fraction of ionizing photons that escapes from 
the halo into the IGM, / csc , and (Hi) the star formation rate (SFR) 
inside the halo. Therefore, we can write: 



^M^ 



1 photons 
M 



M 



s 



(3) 



where the quantity within the brackets is the SFR in a halo of mass 
M, t* is the mean star-formation time-scale in units of the Hubble 
time, iff(z), and /* is the fraction of baryons which is converted 
into stars. For convenience, we combine the free parameters from 
eq. l[3} into a single, redshift-independent, fiducial efficiency of ion- 
izing photons per mass in stars, ead, as: 



x (z) = e fld x 3.8 x 10 s 



H M t H (z) 



photons 

Mq s 



£fid 



0.2 



/csc 

O02 



6.3 x IO^Mq 1 



(4) 



(5) 



These model ingredients are highly uncertain at high redshifts. 
Hence we use this parametrization so that the reader can easily 
substitute his/her own preferred efficiency. The above fiducial val- 
ues were chosen since they provide a reasonably-good fit to high- 
redshift observations: 

1 Note that such extrapolation predicts lower values for the mean free path 
than our choice at z > 9. Thus we likely conservativel y overestimate l m f D 
for the majority of the pertinent redshift range (see also lHaiman et alj200ll) . 
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• The dependence above of the SFR = jfe-M . /*, , 
on the halo mass has been shown to simultaneously re 
produce the obs erved luminosity func t ions of Lyg emittin 



galaxies (LAEs) dshimasaku et al.l 120061; iKashikawa et alj 20061 ; 
IDiikstra et al.l 120071 : Istark et alj|2007l ; iMcOuinn et alj|2007 ), and 
the z = 6 Lyman Break galaxies (LBGs) feouwens et al. 
fairly well, assuming s imilar semi-analytic models (see Fig 7 of 
IDiikstra & Wvifhell2007t) . 

• The fiducial value of / csc =0.02 corresponds to the mean of 
the probability distribution function of / csc that was derived from 
the cumulative Hl-column density distribution along sightlines t o 
(long-duration) 7-ray burst (GRB) afterglows l lChen et alj ;2007). 
T his mean value was also confirmed in the numerical simulations 
of iGnedin etall J2008h . 

• The number of ionizing photons per mass in stars is e 7 = 
6.3 x 10 60 AfQ 1 or equivalently n 7 = e 7 x /im p « 6400 ionizing 
photons per baryon, if o ne assumes a Sa lpeter IMF and a gas metal- 
licity ofZ = O.O4Z dSchaerej|2003h . We point out that for this 
choice of IMF, the number of ionizing photons emitted per baryon 
that is involved in star formation depends on gas metallicity Z as 
n 7 ~ io 4 - 2 - 0029 ^^ 2+9] 2 5 jscnaerer!l2003h . 

Finally, the units of the flux /(x, z) are 'ionizing photons 
s" 1 pcm~ 2 '. On the other hand, the ionizing background is com- 
monly quantified in terms of its mean intensity J(y) which has 
the units 'erg s _1 pcm~ 2 sr _1 Hz -1 '. If the spectrum of the 
ionizing background radiation field is assumed to be of the form 
J{y) = -h\(ylv H Y a x 10~ 21 erg s" 1 Hz" 1 pern -2 sr _1 , then 
the J21 value in the cell located at (x, z) is trivially obtained from 
/(x, z) through the equation 



J 2 i(x,z) = f(x,z) x 



h p a/(An) 



(6) 



10 -21 erg s _1 pem -2 Hz 1 sr _1 
To explicitly show the dependence of our calculated J21 (x, z) 
on the fiducial efficiency choice end, some of our results below will 
be presented in terms of J2i/efid- This allows the reader to sub- 
stitute his/her preferred value of eftd and obtain the corresponding 
physical J21 (x, 2). 

In this study, all calculations assume that only halos with 
masses M > 1.7 x 1O 8 M0 contribute to the generation of the 
ionization and flux fields. This mass at z = 7 roughly corresponds 
to the minimum temperature facilitating efficient atomic cooling, 
10 4 K. However, since this mapping is redshift-dependent, 
this same virial temperature corresponds to somewhat lower halo 
masses at higher redshifts: M = 1.4x 10 8 , 8.6 x 10 7 , and 6.0x 10 7 
Mq at z = 7, 10, and 13, respectively. Because of this, and since 
this atomic cooling cut-off is approximate, we assume a couple of 
different minimum masses when studying the impact of feedback 
on reionization. Note that our fiducial, parametrized UVB distribu- 
tions can be (fairly accurately) adjusted to accommodate different 
contributing minimum masses by adjusting egd accordingly (i.e. 
by the factor of the ratio of the fraction of collapsed matter con- 
tained in halos of mass M > M m j n and our fiducial choice of 
M > 1.7 x 10 s Mq). We stress again that the main goal of this pa- 
per is not to precisely model UV flux distributions at high redshifts 
(such a thing being impossible given our current modest knowl- 
edge/inferences of high- 2 sources); instead we explore a wide pa- 
rameter space to attempt to approximately answer how likely is it 
for radiative feedback to delay the advanced stages of reionization 
when T vir > 10 4 K sources dominate the ionizing photon budget. 
Uncertainties of order unity will thus not deter us from this noble 
task. 



In Figure Q] we present slices through several of our UV flux 
boxes, generated with the above procedure. The slices are shown 
on a linear scale, scaled to the same maximum value. The impact 
of the evolution of structure from z = 13 to 7, and the truncation 
of the flux fields in neutral patches of gas, are qualitatively evident 
in these panels. 

It is important to note that we treat the global neutral fraction 
xhi, redshift z, and the source ionizing efficiency rate per mass 
eion, all as separate free parameters. We can afford this luxury be- 
cause of the speed and versatility of our semi-numerical approach, 
as well as the confidence that our resulting ionization field topolo- 
gies provide a good match t o those generated by cosmolo gical ra- 
diative transfer simulations (Mesinger & Furlanetto 2 0071) . In nu- 
merical simulations of reionization, the source efficiencies are cou- 
pled to the neutral fraction at redshift z. Since we are not modeling 
the progress of reionization in this work (such a task being impos- 
sible at present), we keep ei on and xm separate, the former quan- 
tifying the "instantaneous" ionizing background at z, and the later 
depending on a complicated history of star formation and feedback. 
Again, we do this in order to be as general as possible and explore 
a very wide swath of parameter space. 



3 RESULTS 

3.1 Gas Collapse Fractions 

In Figure [2] we present some results from the spherically symmet- 
ric collapse simulations. The fraction of the gas (normalized to the 
case with no UVBjf] which collapses into halos of mass M at z — 
7, 10, 13 (left to right panels) is shown as a function of the halo 
mass. Curves correspond to J21 = 0.001, 0.01, 0.1, 1.0. The fig- 
ures show that the amount of suppression is a strong function of 
the time elapsed sin ce the halo started being exposed to the UVB 
( IDiikstra et al J2004h . with even 1O 8 M halos at z = 13 being too 
evolved at z on = 14 to lose more ~ 20% of their gas to negative 
feedback even for fluxes as high as J21 = 100 (not shown in the 
Fig.). 

Since our primary interest in this work is how reionization is 
effected by radiative feedback, we will be especially interested in 
the critical specific flux values, JfJ'^M, 2), required to entirelj^ 
suppress the collapse of gas onto halos of mass M at 2. We consider 
this to be more relevant to the issue of reionization than the frac- 
tional suppression of gas seen in Fig. [2] This is because star forma- 
tion likely occurs on a longer time-scale than the time remaining for 
reionization to be completed once a large fraction of the universe 
is ionized and the T v i r > 10 4 halos begin dominating the photon 
budgeQ (se e also for examp le, Fig. 9 in lMesinger et alj|2006l Fig 



budgeO (see also tor example, fig. 9 in Mesinger et al. zUOb, Fig. 
2 in lHaiman & Brvarj|2 006 and Fig. 1 in lLidz et alj|2007l) . The gas 



2 The parameter f gas = M gas (p)/M gas (p = 0), was defined as the mass 
of collapsed gas (i.e. virialized gas) divided by the mass of gas which would 
have collapsed in the absence of pressure. The latter was found to closely 
resemble the J21 = case. 

3 We define "entirely" based on the resolution of our collapse simulations, 
which roughly corresponds to the halo having <0.2% of the gas it would 
have had without being exposed to a UVB. However our results are insensi- 
tive to this particular choice, as the collapsed gas fraction drops extremely 
rapidly with flux in the regime f gas < 10%. 

4 Note that on average star formation should occur on some sizable frac- 
tion (tt, = 2/3 in our model) of the hubble time, tn(z). At Z = 10, 
(2/3)tj/ = 500 Myr; this time-scale is larger than, for example, the time 
elapsed from 2 = 10 to 2 = 7, At = 300 Myr. 
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Figure 1. Slices through our UV flux boxes shown on a linear scale, scaled to the same maximum value. Slices correspond to (z, £hi) = (7, 0), (10, 0), (10, 
0.63), (13, 0). Each slice is 0.5 Mpc thick and 100 Mpc on a side. 




M/M e M/M s M/M Q 



Figure 2. Fraction of gas (normalized to the J21 = case) which collapses into halos of mass M at 2 = 7, 10, 13 (left to right panels). Curves correspond 
to J21 = 0.001, 0.01, 0.1, 1.0, top to bottom. All runs assume 2on — 14. 



which was prevented from collapsing in Fig. [2] was in the outer 
regions of the halo, and would likely have formed stars later, prob- 
ably after reionization was completed. Halos which retain even a 
little of their gas reservoir can form stars fairly quickly, as this gas 
is close to the high density core where cooling is efficient. Thus 
the contribution of halos to the "relative immediacy" of complet- 
ing reionization is better judged with J%1 (M, z). However, since 
these assumptions are by no means certain, we present total feed- 
back estimates in i)3.3l based on two models: (i) partial gas suppres- 
sion, in which we use the collapsed fraction, /gas, that is plotted in 
Fig. [2] and (ii) and total gas suppression, in which halos that see a 
UVB that is greater (less) than J cr i t retain none (all) of their gas. 

In Figure[3]we show the critical Jf™ (M, z) values required to 
entirely suppress the collapse of gas onto halos of mass M at z = 7 
(bottom curve) and z = 10 (top curve). We obtain Jff (M, z) 
at each M by extrapolating the results from our collapse simula- 
tions linearly in / gas vs. ln(J2i). It is clear from the figure that 
J21 (M, z) is a strong function of halo mass. At this stage, one 
can already estimate that the contribution of ~ 10 s Mq halos to 
reionization could be suppressed in the later stages, but suppress- 
ing > few x 10 8 Mq halos would be extremely difficult, even given 
our conservative assumptions. 



3.2 UV Flux Distributions 

Following the procedure outlined in i]2.2.1| we generate UV flux 
fields at z = 7, 10, and 13, and for several values of Shi. Proba- 
bility distributions of J21 /efld extracted from these semi-numerical 




1O-4 1 1 1 1 — 11111 

10 8 10 9 

M/M 



Figure 3. J21 values required to entirely suppress the collapse of gas onto 
halos of mass M at z = 7 (bottom curve) and z = 10 (top curve). 

simulations are displayed in Figure [4] Dotted curves are generated 
from all pixels whereas solid curves are generated only from pixels 
which contain halos. The dashed curves in the z = 10 panel are 
generated from pixels which contain halos, assuming xm ~ 0, but 
with a different m.f.p., l m f p = 30 Mpc. 

Comparing corresponding dotted and solid curves clearly 
shows that halos are irradiated on average by an enhanced ioniz- 
ing background, due to their clustering. This enhancement is larger 
at early times when sources are more biased. Thus at z — 13 the 
peaks of the distributions generated from halo locations are a fac- 
tor of ~4 larger than the peaks of the distributions generated from 
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z=7 z=10 z=13 




Figure 4. Probability density functions (PDFs; top) and cumulative probability density functions (CPDFs; bottom) of J21 Afid at z = 7, 10, 13 (left to right). 
Dotted curves are generated from all pixels whereas solid curves are generated only from pixels which contain halos. Curves correspond to xm ~ 0, 0.31, 
0.64, 0.81, 0.94 (right to left). The dashed curves in the 2 = 10 panel are generated from pixels which contain halos, assuming xm ~ 0, but with a different 
m.f.p., i m f p = 30 Mpc. 



random locations in space, with this factor decreasing to <1.5 by 
z = 7. One can also notice that as sources become less biased at 
lower redshifts, the flux distributions at low xm become narrower 
and more gaussian. This is due to the fact that a larger fraction of 
the newly formed halos are born in less-biased regions with a corre- 
spondingly weaker UVB. These new arrivals thus push the low-flux 
tail to higher values and narrow the distributions. 

The imprint of our choice for the m.f.p., l m f p , is also seen in 
the high-flux tails of the distributions in Fig. [4] especially in the 
solid curves. There are two regimes governing the high-flux tail: 
(1) at low xhi < 0.5, the m.f.p sets the maximum flux seen by ha- 
los; (2) at high Shi > 0.5, the HII bubble size sets the maximum 
flux seen by halos. Thus the high-flux tails of the distributions over- 
lap at low xhi < 0.5 [regime (1)]; then at im > 0.5 the tails of 
the distributions drop and gradually decrease as the bubble size de- 
creases and the halos see fewer and fewer sources [regime(2)]. The 
transition at xm ~ 0.5 denotes the moment when the characteristic 
bubble size (specifically the large tail of the distribution where the 
clustered sources lie) surpasses the Z m f p ; this is confirmed in Fig. 
[5] where we plot the volume-weighted size distributions of ionized 
region^ 

Also shown in the middle panel of Fig.|4]are the flux probabil- 
ity distributions (clashed curves) created from pixels which contain 
halos, assuming xm ~ 0, but with a 50% higher m.f.p., Z m f p = 30 
Mpc. By comparing these curves to the right-most solid curves, we 
can isolate the maximum impact of Z m f p , when bubble size is not an 
issue, i.e. xm ~ 0. This increase in the m.f.p. shifts the peak of the 
flux distribution by ~ 30%; however, the shift is less pronounced 
at the high-flux than the low-flux end, since the ionizing radiation 
from nearby halos dominates over the large-scale background at the 
high-flux tail. This, combined with the fact that the high-flux, high 

5 These distributions were created by randomly choosing an ionized cell, 
and recording the distance from that cell to a neutral cell along a randomly 
chosen direction. We repeat this Monte Carlo procedure 10 7 times for each 
Shi- Note that the ionization topology and bubble d istributions at fixed 
xm are not sensitive to redsh ift over this interval dMcOuinn et alj[2007l ; 
iMesinger & F urlanetto 2008b). 




R (Mpc) 

Figure 5. Volume-weighted size distributions of ionized regions at z = 10 
and x m us 0.94, 0.81, 0.63, 0.45, 0.31, 0.11 (left to right). The vertical 
dashed line denotes our fiducial choice of i m f p = 20 Mpc. 

xm tails overlap more at high z, when sources are more biased, 
suggests that the high-flux tails can effectively be modeled by a 
small-scale, Poissonian halo clustering component combined with 
a smooth contribution from larger scales (c.f. Dijkstra et al. 2008, 
in preparation). Similarly, the high-flux tail of the UVB distribu- 
tions might be extended if stocasticity in the source efficiencies is 
taken into account (Dijkstra et al. 2008, in preparation); however, 
in this paper, we study the impact on reionization, and any signifi- 
cant impact would have to result from fluxes not too far out on the 
high-value tail. 



3.3 UV Radiative Feedback 

Combining our parametrized flux boxes with the results from our 
collapse simulation, we can proceed to answer the fundamental is- 
sue of this paper: can UV radiative feedback substantially suppress 
the contribution of T v i r > 10 4 K halos to the advanced stages of 
reionization, when such halos are expected to dominate the ioniz- 
ing photon budget? We do this by counting the fraction of halos 
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Figure 6. Top: Fraction of M > 1.7 X 10 Mq halos left without gas at z = 7 and 10 {left to right). Bottom: Same as the top panels, but weighted by each 
halo's total mass. Curves correspond to egd = 0.25, 0.75, 2.25, 6.75, 60.8 {bottom to top in the left panels) and egd = 25, 75, 225, 675 {bottom to top in the 
right panels). 



of mass A/, which lie in regions of flux stronger than the critical 
flux, Jf" {M, z), shown in Fig. [3] assuming various choices for 
our ionizing efficiency parameter, eea- 

In the top panels of Fig. [6] we show the fraction of M > 
1.7 x 1O 8 M halos left without gas at z = 7 and 10 (left to right). 
The bottom panels show the same quantity, but weighted by each 
halo's mass. These mass-weighted gas fractions are directly pro- 
portional to total gas mass 'locked up' in dark matter halos, which 
in turn determines the global star formation rate. The bottom pan- 
els therefore present a more realistic picture of the impact of UV 
feedback than the top panels. 

From Fig. [6] it is evident that "reasonable" choices of ionizing 
efficiencies, efld ~ 1 are insufficient to quench star-formation in 
the majority of halos. At z = 7 (10) efficiencies of end ~ 5 (200) 
are required to completely evacuate/suppress gas from roughly half 
of M > 1.7 x 10 s halos. The required efficiencies become 
even more out-of-reach when considering the mass-weighted gas 
fractions. For example, suppressing more than half of the global 
gas reservoir requires egd ~ 100 (1000) at z — 7 (10)! In our 
scenario, halos at z = 7 are much more susceptible to radiative 
feedback than those at z = 10, since they have been exposed to an 
ionizing background radiation for twice as long (but note that we 
have exaggerated this effect, see § 12. It , 

We also note that the feedback effects are not particularly sen- 
sitive to the ionization efficiency, egd, once that parameter goes be- 
yond a certain minimum value (e.g. efld > 1 at z — 7 or efld > 100 
at z — 10). This can be understood by noting that the critical flux 
curves in Fig.[3]are very strong functions of mass. Once ionizing ef- 
ficiencies are large enough that our smallest mass scale experiences 
feedback, additional large increases in efld are required to signifi- 
cantly suppress gas collapse onto more and more massive halos. 

One can note the footprint of i m f p in Fig. [6] as well. The rise 
in the curves becomes shallower at the lower neutral fractions, 



xhi < 0.5, when the characteristic bubble size surpasses l m t p , as 
discussed above. 

It is useful to remove the mass function from the cumulative 
feedback effects depicted in Fig. [6] and look at radiative feedback 
as a function of mass scale. Thus, in Fig. [7] we plot the fraction of 
halos of mass M left without gas at z — 7 and 10 (bottom to top). 
Solid curves correspond to xm. ~ 0, while dotted curves corre- 
spond to Shi ~ 0.81. We see that roughly an order of magnitude 
increase in the ionizing efficiency is required to extend the mass 
scale susceptible to strong feedback by a mere factor of ~ 2. Only 
M < 2 x lO 8 A/0 halos at z = 7 experience significant feedback 
effects at our fiducial efficiency efld ~ 1. 

As discussed above, it is likely that estimates based on 
J|" 4 (M, z) are decent indicators of the strength of feedback ef- 
fects during the advanced stages of reionization. We justified this 
hypothesis in ^3. 1 1 with a crude comparison of the probable time- 
scales of star formation and the completion of reionization. This 
"all or nothing" scenario assumes that all halos of mass M exposed 
to specific fluxes in excess of J%1 (M, z) no longer contribute ion- 
izing photons, whereas those exposed to weaker fluxes continue 
doing so unhindered until the end of reionization. Nevertheless, it 
is useful to also consider a model with a less sharp cut-off due to 
feedback. In this extreme case, we assume that the halos are able 
to "rapidly" (i.e. before the end of reionization) convert their entire 
gas reservoir at z to ionizing radiation. In this case, radiative feed- 
back should follow the relative suppression of collapsed gas shown 
in Fig. [2] Additionally, from a more fundamental standpoint, it is 
interesting to note the extent to which a UVB depletes the global 
gas reservoir at high-jz. 

Therefore, in Figure[8]we plot the fraction of gas (normalized 
to the J21 = case) still remaining in M > 1.7 x 10 s M Q (solid 
curves) and M > 10 s Mq (dashed curves) halos at 2 = 7 and 10 
(bottom to top). Pairs of curves correspond to egd = 0.1, 1, 10, 100 
(top to bottom), at both redshifts. 
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Figure 7. Fraction of halos of mass M left without gas at z = 1 and 10 
(bottom to top). Solid curves correspond to xjji ~ 0, while dotted curves 
correspond to ipjl ~ 0.81. Each pair of curves assumes egj = 250, 25, 
2.5 (right to left in the top panel), and Efld = 25, 2.5, 0.25 (right to left in 
the bottom panel). 

As in the top curves of Fig. [6] feedback is not very sensitive 
to esd, requiring large increases in efld to suppress gas in more and 
more massive halos. This is especially evident at z = 7, where 
feedback effects are strongest: here we see the remaining gas frac- 
tion only evolves < 20% as efld is increased from 0.1 to 100. An 
effective "asymptote" at / gas ~ 0.4 is reached as the majority 
(though not all) of the gas is prevented from collapsing onto the 
"vulnerable", 10 8 < M < 10 9 halos. The suppressed fraction, 
1 — /gas ~ 0.6, roughly corresponds to the global collapse frac- 
tion of these halos at z = 7. This effective asymptote is at a lower 
/ gaa at z = 10, reflecting the smaller characteristic halo mass at 
that redshift and the corresponding higher global collapse fraction 
contained in these small, "vulnerable" halos. 



4 ASSUMPTIONS AND UNCERTAINTIES 

In ^ 12. II we investigated the impact of a photoionizing flux on the 
ability of the gas to cool and collapse onto dark matter halos by per- 
forming 1-D hydrodynamical simulations. These calculations over- 
estimated the impact of a UVB on the ability of gas to cool and 
collapse into the centers dark matter halos for a number of reasons: 

• We assumed the spectrum of the UVB to be hard: J(v) oc 
J21 {v/vh) -01 and a = 1.0, which is typically associated with 
accreting black holes. However, stellar sources are likely to sig- 
nificantly contribute to the UVB at these r edshifts. These sources 
have much softer spectra, with a ~ 5 (e.g. lBarkana & Loeb 2001, 
p.70). Under the assumption that the gas is transparent to ion- 
izing photons, the photoionization [photoheating] rate scales as 
oc J21/O + 3.0) [cx J 2 i/(q + 3.0)(q + 2.0)]. Choosing a hard 
spectrum is (almost) equivalent to boosting J21 , because in our sim- 
ulations self-shielding is ignored and the gas is transparent to ion- 



Figure 8. Fraction of gas (normalized to the J21 = case) still remaining 
in M > 1.7 X 10 8 M Q (solid curves) and M > 10 8 M Q (dashed curves) 
halos at z = 7 and 10 (bottom to top). Pairs of curves correspond to egd = 
0.1, 1, 10, 100 (top to bottom). 

izing photons. Thus if the UVB is dominated by stellar sources at 
these high redshifts, we expect UVB intensities to be a factor of ~ 
five lower than our conservative estimates (effectively end ~ 0.2 
for our reionization estimates)! 

• Our simulations ignore that at sufficiently large densities (5 > 
100, e.g. Yang et al, 2006), the gas is capabl e of self-shielding 
again st the ionizing background (also see e.g. ISusa & Umemural 
2004). The reasons for ignoring self-shielding in our calculations 
are ((') it speeds up the simulations: properly accounting for self- 
shielding requires one to compute the increasing hardness of the 
spectrum as one moves deeper into the cloud, and (») in more re- 
alistic 3-D models of collapsing gas clouds, ionizing photons may 
penetrate deeper into the cloud via paths of lower HI column den- 
sities than one expects from a 1-D calculation. Hence, if one ac- 
counts for self-shielding using a 1-D simulation, one may under- 
estimate the impact of the UVB, while throughout this paper we 
have conservatively chosen to overestimate the impact of the UVB. 
For example, Dijkstra et al (2004) found that ~ 95% of the gas 
was able to collapse into dark matter halos with u c i rc = 15 km s -1 
(roughly corresponding to 10 8 A/q) that collapsed at z = 11 for 
Zon = 17 and J21 = 10~ 2 when self-shielding was included. This 
is a significant rise compared to the ~10% fraction that was found 
for a model in which self-shielding was ignored. 

A large uncertainty in our determinations of the UVB is the 
ionizing photon mean free path inside ionized regions, Z m f p . Our 
choice of Z m f P = 20 Mpc is consis t ent with other theoretic al es- 
timates (e.g. iFurlanetto & Ohll2005t IChoudhurv et alj|2008h . and 
might even be conservatively large (i.e. lowering l m f p only serves 
to reduce the calculated values of the UVB) at 2 > 9 given the 
extrapolation of LLS density f rom lower redshift ob servations: 
Z mfp w 110[5/(1 + z)] 3 Mpc jMiralda-Escudl 120031) . However, 
given that that the nature of these absorbers is unknown, all cur- 
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rent estimates are just educated guesses so similarity among them 
is a somewhat underwhelming validation. Nevertheless, we demon- 
strated that Z m f p only affects the UVB when the characteristic bub- 
ble size (weighted by the halo occupancy) surpasses i mlp . There- 
fore, if Z m f p is larger, then we have underestimated the value of the 
UVB, but only during the very late stages of reionization, when the 
characteristic bubble size exceeds l m t p > 20 Mpc. Hence, we con- 
clude that our main conclusion that UV feedback does not delay 
the bulk of reionization is completely unaffected by the uncertainty 
in Z m f p . 

Although high-redshift sources remain poorly understood, our 
fiducial choice of source ionization efficiency, e^, was moti- 
vated by its ability to r eproduce the observed luminosity func- 
tions of z ~ 6 LAEs ([Shimasaku et al.l 120061 : iKashikawa et al.l 
|2006|) . and LBGs dBouwens et alj2006l) fairly well, assuming sim- 
ilar semi-analytic m odels (Diik stra et al J 120071 : IStark et alj 2007; 
iMcO uinn ct al ]|2007l) . It seems that significantly delaying the bulk 
of reionization through UV feedback is only likely if ionizing ef- 
ficiencies of z > 10 sources are much higher (a couple of orders 
of magnitude; see right panel of Fig.[6j than such z ~ 6 data seem 
to imply. An factor of ~ 4 can be squeezed-out from eq. 10 if one 
assumes all star formation occurs on a dynamical time, and an addi- 
tional factor of ~ 10 can be obtained if one a ssumes ionizing emis- 
sivity, 6y, expected of Population III stars (Schaerer 2003). This 
seems unlikely, although such a strong redshift e volution in ioniz- 
ing efficiency is implied if the Stark et al. (2007) lensed sample of 
z ~ 9 LAEs is genuine dStark et al.l l2007; Mesi nger & Furlanettd 
2008b). Furthermore, efld depends linearly on / csc , which we as- 
sumed to be / csc = 0.02 (see § 12.2. It . To our knowledge, no ob- 
servation s exist that hint at significantly larger average values o f 
/csc (e.g. llnoue et alfcoodlshaplev et alj|2006l : lsiana et alj|2007h . 
Theoretically, escape fractions of order unity (/ osc ~ 1) only ap- 
pear associated with minihalos that contain a single massive star 
i lWhalen et ai1l2004l) . 

Again, we point out that we do not model reionization and 
feedback self-consistently. Our conclusions on the strength of feed- 
back effects are thus conservative, in that results from WMAP sug- 
gest most halos would have been exposed to a UVB for shorter 
periods of time than our choice of z on = 14 implies (see discus- 
sion in § I2.lt . We are also conservative in neglecting the build- 
up of the UVB from z on — 14, assuming a constant distribution 
calculated at z throughout the entire interval z on — > z. On the 
other hand, as we point out above, we assume a minimum mass 
of M m i n > 1.7 x 10 s Mq when generating ionization and flux 
fields, roughly corresponding to T v i r = 10 4 K at z = 7. But the 
same virial temperature corresponds to a somewhat lower M m i n 
at higher redshifts. If star formation is indeed allowed down to 
exactly T v i r = 10 4 K, the contribution of neglected halos with 
M < 1.7 x 1O 8 M0 can be approximated with a ~ 40% increase 
in the effective efficiency egd at z = 10 (assuming the ionizing 
photon contribution traces the collapsed fraction, as was done in 
this work). This uncertainty is smaller the other uncertainties in 

Efld- 

Even assuming our source ionizing luminosity prescriptions 
are realistic, we do not expect to accurately model the tails of the 
UVB distributions, due to the following effects: 

• Cosmic variance on the scale of our simulation box, 100 Mpc, 
causes us to miss statistically rare regions. 

• Our halo fields are generated on a 1400 3 grid and our velocity 
fields are generated on a 700 3 grid. Thus we cannot model the halo 
locations to better than 100 Mpc/700 » 0.14 Mpc. This means that 



the very rare, very close source configurations are missed by our 
code. 

• We don't include a "duty cycle" in our source prescriptions, 
assuming that all sources are "turned-on". While our prescription 
in eq. l|3} accounts for this effect on average by sweeping it under 
the effective /» or / CS c parameters, the duty cycle stocasticity can 
extend the tails of the UVB distributions (Dijkstra et al. 2008, in 
preparation). A similar effect is also produced if the absorbers are 
treated as distinct objects, instead of merely accounting for their 
mean flux decrement, exp(— r/£ m f P ), as in eq. (fTJ. 

Nevertheless, the tails of the UVB distribution are not impor- 
tant in answering the pertinent question of "does UV radiative feed- 
back significantly affect the progress of the bulk of reionization". 
Ionizing fluxes close to the means of the UVB distributions are 
more relevant in answering this question. Given our generalized 
parametrizations, we are confident our models are sufficiently ac- 
curate for these purposes. Thus we do not expect to loose sleep 
worrying about the effects mentioned in the previous paragraph. 



5 CONCLUSIONS 

In this paper, we quantify the importance of UV radiative feed- 
back during the middle and late stages of reionization. Specifi- 
cally, we concern ourselves with halos capable of atomic cooling, 
T v i r > 10 4 K. We first run suites of spherically-symmetric halo 
collapse simulations using various values of z, M, and J21. We 
then generate parametrized UV flux distributions at z = 7, 10, and 
13, using semi-numerical, large-scale simulations of halo and ion- 
ization fields. Combining these two results, we estimate how effi- 
cient is radiative feedback at hindering the progress of reionization, 
during its advanced stages (when T v i r < 10 4 K halos become sub- 
dominant contributors of ionizing photons). This tiered approach 
allows us to explore an extremely wide range of parameter space, 
which is necessary to make any robust conclusions, given our poor 
knowledge of the properties of high- z sources. 

We find that under all reasonably conservative scenarios, UV 
feedback on atomically-cooled halos is not strong enough to no- 
tably delay the bulk of reionization. For fiducial choices of source 
ionizing efficiencies (calibrated to match z ~ 6 LAE and LBG 
luminosity functions) and turn-on redshift, fewer than 40% of 
T vir > 10 4 K halos are left without gas at (z, xm) ~ (7, 0); this 
number drops to ~ 10% when the distribution is mass-weighted (as 
is more appropriate for estimating the global star formation rate). 
Suppressing more than half of such halos requires a factor of ~ 
3^1 increase in fiducial ionizing efficiencies at (z, xm) ~ (7, 0) 
and over two orders of magnitude increase for the same fraction at 
{z,x m ) «(10,0)! 

We find that feedback is very strongly dependent on halo mass. 
For example, at (z, xm) ~ (7, 0) and efld = 0.25, the fraction of 
halos left without gas decreases from 1 to 0.2 as the halo mass is 
increased only from 10 8 to 1.7 x 10 8 Mq. The fraction of affected 
halos only decreases as xm is increased. 

Accurate quantitative estimates will have to wait for a break- 
through in our understanding of the UV emission properties of 
high-redshift sources, and their dependence on halo mass. Never- 
theless, the inability of photoionization feedback to delay the mid- 
dle and late stages of reionization is compelling, especially given 
that we do not include radiative transfer effects which would only 
decrease the relevance of UV feedback (see §|4]for a more detailed 
discussion). 

It seems that delaying the advanced stages of reionization 
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through UV feedback on M > 10 s Mq halos is only likely if ion- 
izing efficiencies of z > 10 sources are much higher (~ two orders 
of magnitude) than z ~ 6 data seem to imply. An evolution in end 
could be obtained by changing the fiducial values in eq. {4j, possi- 
bly as a result of a top-heavy IMF or a star-burst dominated epoch. 
However, an increase in the fiducial efficiency, end, by ~ two or- 
ders of magnitude seems unlikely. 

Our results also suggest that the natural time-scale for the bulk 
of reionization is the growth of the collapsed fraction contained in 
T v i r > 10 4 K halos. The natural timeframe for half of the uni- 
verse to be ionized in such a prescripti on, using a toy model for 
the evolution of the HII filling facto r (Furlanetto & Loeb 2005; 
IMesinger et all 120061; lHaiman & Brvanll2006l : llidz et al.ll2007l) is 
Azre ~ 2 - 3, with a small filling factor tail extending to higher 
redshifts. This tail is governed by more complicated feedback on 
T vir < 10 4 K objects. Such a scenario would be consistent with 
results implied by both the WMAP optic al depth measurements 



(Dunklev et aHl2008l: 



(Mesinger & Haiman 



iMesinger & Haimanl 



Komatsu et al. 2008) and the SDSS quasars 
2004IWvithe & Loebll2004 iFan et alj|2006l : 



20071) . especially given the inhomogeneous 
and the associated difficulties in inter- 



nature of reionization 
preting present data dlidz et al.l 120061; iBolton & Haehneltl |2007|; 
Masel li et alj 120071 ; iBecker et alj |2007| ; IMesinger & Furlanetto! 



2008al) . This scenario is also in agreement with recent observational 



claims that faint galaxies at or below current de tection thresholds 
domin at e the ionizing photon budget at z > 6 {Yan & Windhorst 
l2004al lbl; lBouwens et all2006h . 

Finally, our results underscore the importance of extended dy- 
namical ranges when modeling reionization. Simulations must be 
capable of resolving halos with mass > 1O 8 M0, even when mod- 
eling the late stages of reionization, while at the same time being 
large enough to capture HII regions several tens of megaparsecs in 
size. 



We are grateful to Steven Furlanetto and Zoltan Haiman for 
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